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Abstract 


This study investigates the function of 5-HT type 3 (5-HT3) receptors using a computational approach. Antagonists of the 
5-HT3 receptor are currently one of the most effective therapeutic agents in treatment of chemotherapy-induced nausea, 
vomiting, and irritable bowel syndrome. Several experimental studies have shown the effect of pharmacological agents 
such as 9-tetrahydrocannabinol (THC), the psychoactive component of Cannabis, Cannabidiol (CBD), a non-psychoactive 
ingredient of Cannabis plant, Menthol, Propofol, and etc. on the functional human 5-HT3 receptors expressed in Xenopus 
oocyte as well as rat nodose ganglion neurons. 5-HT evoked currents recorded by a two-electrode clamp technique were 
inhibited by ligands in a concentration dependent manner. Simulations of allosteric inhibition were modeled using Vina 
docking techniques with the 5-HT3 structure. The 5-HT3 structure was found using homology sequence similarity 
techniques with the neuronal nicotinic acetylcholine receptor (nACH) and inhibitory neurotransmitter receptor for GABA 
(A). Results of studies with other members of the super family of ligand gated ion channels signified key residues 
involved in ligand binding sites within the transmembrane region of 5-HT3. Flexible and rigid docking simulations 
around key residues resulted in a number of low-energy (high affinity) configurations of ligand binding. The predicted 


residues TYR and THR may constitute a naturally occurring binding site for 5-HT3. 


Introduction effective therapeutic agents in treatment of chemotherapy 


induced nausea, vomiting, and irritable bowel syndrome 
[1]. The serotonin type three (5-HT3) receptor, a member 
of the cys-loop family of ligand-gated ion channels, 


Plants have long been used as a source of natural 
chemicals for medicinal products on a commercial scale. 


Menthol, a chemical constituent of the mint plant is used 
as a very important flavoring agent in products such as 
toothpaste, mouthwash, foods and cigarettes. Studies 
have shown that menthol containing products can be used 
for medicinal purposes ranging from pain management to 
nausea and vomiting treatment. Until recently, the cellular 
and molecular targets mediating beneficial effect of 
menthol containing products was unknown. Antagonists 
of the 5-HT3 receptor are currently one of the most 
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mediates fast depolarizing actions of 5-HT in the central 
and peripheral nervous system [2]. 


Several experimental studies have shown the effect of 
pharmacological agents inhibit the membrane currents of 
5-HT. The effects of 9-tetrahydrocannabinol (THC), the 
psychoactive component of cannabis, on the function of 5- 
HT type 3 (5-HT3) receptors were investigated using a 
two-electrode voltage clamp technique in Xenopus 
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oocytes, and a whole-cell patch clamp technique in rat 
nodose ganglion neurons [3]. The results of this study 
indicate for the first time that, in Xenopus oocyte 
expression system, the inhibition by THC, a major 
psychoactive constituent of Cannabis sativa, of 5-HT3 
receptors depends significantly on the amount of 
functional receptors expressed in this system. Also, the 
same laboratory investigated the effects of alcohol 
monoterpene menthol, a major active ingredient of the 
peppermint plant, on the function of human 
hydroxytryptamine type 3 (5-HT3) receptors expressed 
in Xenopus oocytes [4]. They found that 5-HT3 receptor 
mediated currents in acutely dissociated nodose ganglion 
neurons were also inhibited by menthol (100 uM). These 
data demonstrate that menthol, at pharmacologically 
relevant concentrations, is an allosteric inhibitor of 5-HT3 
receptors. 


However, experimental studies don’t include the 5-HT3 
receptor kinetics such as binding sites as well as binding 
energy in their electrophysiology data. Clearly, more 
detailed investigations in mammalian neurons or cell 
lines are needed, where the receptor kinetics in a more 
detailed and precise manner are included. In this 
computational study we try to investigate the function of 
5-HT type 3 (5-HT3) receptors using same 
pharmaceutical agents in experimental data [5]. We 
simulate this experiment using docking of small ligands 
onto the 5-HT3 chain a receptor (Figure 1). 
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Figure 1: 5-HT3 chain a receptor. 


Docking predicts the position and orientation of a 
ligand when it is bound to a protein receptor or enzyme. A 
scoring function is used to calculate the energy of 
interaction based on electrostatics, Van der Waal’s 
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interactions, hydrogen bonds, solvation effects, and loss of 
entropy. 


lon pore 


Methods 


Experimental Methodology 


Synthesis of cRNA from cDNA clones of the human 5- 
HT3A, Radioligand Binding Studies, [35S] GTPgS Binding 
and Whole-Cell Patch-Clamp Recording were performed 
[4,6]. 


Computational methodology 


e Obtain a Structural Model of 5-HTR3A. We used 
Homology modeling techniques utilizing the PHYRE2 
6 server to generate a structural model for 5-HTR3A. 
The result was a model for each chain of the 
pentameric structure of 5-HT3 [7]. Chain A of the 
nicotinic acetylcholine receptor (nAChR), under PDB 
code 2BG9 and the g-aminobutyric acid a receptor 
(GABAAR), PDB code 4COF was chosen as the 
modeling template for a Cys-loop receptor. The 
transmembrane domain of the modeled structure was 
employed for docking, as described below. 

e Docking of Ligands on the 5-HTR3A Structural Model. 
We use Autodock Vina program to estimate binding 
energies and binding modes of ligands such as 
menthol, eucalyptol, propofol, cyclohexane, THC, and 
CBD. *** and thymol and 2isopropylphenol*** In 
order to use this docking program we have to define 
the region of the receptor related to a binding site. 
This region is identified based on experimental 
studies1,2,3. Autodock Vina performed rigid and 
flexible docking simulations. 


Choice of the specified grid box area took into account 
similar binding characteristics believed to be shared by 
propofol and menthol [8] and the close homology of the g- 
aminobutyric acid a receptor (GABAAR) to 5-HTR3A [1]. 
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Flexible docking. A more sophisticated model of 
receptor allows close residues inside of binding site to 
adjust side chain conformations. This was performed on 
all large ligands (THC, CBD) 

Open macromolecule file with Harlem and assign binding 
site. 


Find center of bind site and measure a sphere 15A 
away from center. Color and note each amino acid that 
lies within the measured area. Assign amino acid residues 
within site for flexible binding techniques in autodock. 


e Analyze Interactions. Each ligand was modeled and 
interactions were visualized using Autodock Tools. 


First Open ‘Autodock tools’ > click ‘AutoDock vina 
Result’ 

Locate receptor and add as macromolecule + select ‘vina 
result’ 

Macromolecule - vina_ result 

(ligand) - Propofol.pdbqt 

Click ‘Analyze Results’ and select ‘Single molecule with 
multiple conformations’ click ‘OK.’ 


5HT3.pdbqt, 


The docked molecule has now displayed its lowest 
energy conformation. Keyboard left and right to browse 
other vina docking conformations lowest to highest. (9 
total - all zoomed out w/out interactions with amino acid 
chain of 5HT3) 


Now click ‘Analyze’ > ‘Dockings’ > ‘Show Interactions’ 
(This changes the display to a close up image of the first 
lowest conformation of the docked ligand.) 


The background color was set to white. (Customizable 
within display widget) 


Click ‘Show Molecular Surface’ for ligand and color by 
atom type to find ligand atoms close to receptor atoms 
and white for atoms not close to any receptor atoms. Click 
‘Create Spheres’ to show wire spheres around receptor 
atoms close to a ligand. (color gray) 


Add labels and sticks for residues interacting close to 
ligand. (optional: choose within 5-10A) 


Click ‘Hydrogen Bonding’ - choose ‘Activate as Spheres’ 
(color green - default) and adjust radii + spacing for 
optimum view. 

Toggle and zoom to get preferred angle that clearly 
captures ligand atom bonding with amino acid residues or 
hydrogen bonds, or both. Capture image by clicking ‘Set 
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Antialiasing’ which will increase quality of viewing 
window and click save image. 


Output list of ligand contacts, receptor contacts and 
hydrogen bonds. 


Results 
Experimental Results 


Comparison of the effects of equimolar concentrations 
(100 uM) of cyclohexane, eucalyptol, propofol, and 
menthol on 5-HT3 receptors. Bars represent the means + 
S.E.M (n=5-8). Chemical structures of the compounds 
were presented at the top of each bar [4]. 
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Computational Results 


A. Results of rigid docking to lower binding site next to 
Threonine 361 (menthol binding site, lowest energy 
conformation) [7]. 


Compound Binding energy (kcal/mol) 
Propofol -6 
Menthol -5.6 
Eucalyptol -4.5 
Cyclohexane -4.5 
Compound Binding energy (kcal/mol) 
Thymol -6.4 
2-isopropylphenol -6.1 
Benzene -6 
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Compound 2-D Ligand interaction diagram 


3-D depiction of interaction 


Propofol 


Menthol 
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Eucalyptol 


Cyclohexane 


B. Results of rigid and flexible docking to upper binding The following are readout files for the rigid docking of 
site near Tyrosine 346 (THC and CBD binding site, THC and CBD. The docked molecule has now displayed its 
lowest energy conformation). lowest energy conformation. 
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Hae aE aE aT aE aE aE ae ae aE aE aE aE aE aaa aaa aaa aaa aaa 


# If you used AutoDock Vina in your work, please cite: 
O. Trott, A. J. Olson, 

with a new scoring function, efficient optimization 
455-461 


DOI 10.1002/jcc.21334 


Output will be €hé 6ut.pdbaqt 


AutoDock Vina: improving the speed and accuracy of docking 


multithreading, Journal of Computational Chemistry 31 (2010) 


Please see http://vina.scripps.edu for more information. 
HEHEHE HH HH HH HH EE EE EE EE EEE EE HE HH EE EE EE EE EE EE EEE EE EE EE EE EE EE EE EE 


and 
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Detected 4 CPUs 

Reading input ... done. 

Setting up the scoring function ... done. 

Analyzing the binding site ... done. 

Using random seed: 1209244732 

Performing search ... done. 

Refining results ... done. 

mode | affinity | dist from best mode 
| (kcal/mol) | rmsd 1.b.| rmsd u.b. 
----- +------------4----------4+---- 
1 8.4 0.000 0.000 
2 -8.4 1.768 6.373 
3 —8:.3 994 2.499 
4 -8.2 1.816 6.953 
5 -8.0 1.698 2.440 
6 =e 958 7.144 
7 Ed 1.816 6.167 
8 -7.6 960 6.766 
9 -7.6 1.946 2.376 
Writing output ... done. 


HEHEHE HH HEH HH EH EE EE EE EE HEH EE HE EH EE EE EE EE EE EE EEE EE EE EE EE HE EE EE EE 
If you used AutoDock Vina in your work, please cite: 


O. Trott, A. J. Olson, 

AutoDock Vina: improving the speed and accuracy of docking 
with a new scoring function, efficient optimization and 
multithreading, Journal of Computational Chemistry 31 (2010) 
455-461 


DOI 10.1002/jcc.21334 


Please see http://vina.scripps.edu for more information. 


Hae aE aE aT aE aE aE aE ae aE aE aT aE aE ae aa aaa aaa aaa aaa 
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Output will be ebd 6ut.pdbaqt 


Detected 4 CPUs 

Reading input done. 

Setting up the scoring function done. 

Analyzing the binding site done. 

Using random seed: 1410998816 

Performing search done. 

Refining results done. 

mode | affinity | dist from best mode 
| (kcal/mol) | rmsd 1.b.| rmsd u.b. 
----- +------------4----------+4---------- 
1 --7.7, 0,000 0.000 
2 = Tb 0.223 1.863 
3 -7.2 2.642 Sel oy 
4 -7.2 2.51.0 5.303 
5 -6.2 6.2175 12.965 
6 -6.1 1.014 2.944 
7 -6.0 2.121 5.743 
8 -6.0 25S 5.161 
9 =535 5.423 11.766 
Writing output done. 


C. Flexible Docking 


Flexible docking 


between 


Delta-9- 


Tetrahydrocannabinol (THC), Cannabidiol (CBD) and 
residue receptor (5-HT3) - flexible residues include TYR 
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251, LEU 255, LEU 301, THR 304, VAL 345, CYS 350 and 


LEU 509. 


Rigid and flexible docking results for THC and CBD and 


their corresponding 
displayed in 3D below. 


lowest energy conformation 


Compound Rigid docking - Binding energy | Flexible docking - Binding Energy 
(kcal/mol) (kcal/mol) 

Delta-9-Tetrahydrocannabinol -8.4 -9.2 

Cannabidiol -7.0 -8.8 


THC = -9.2 kcal/mol 


CBD = -8.8 kcal/mol 
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Pi-Pi bonding interactions between Delta-9-Tetrahydrocannabinol (Ligand) and residue PHE 347 receptor (5-HT3) - 
attractive, noncovalent interactions between aromatic rings. (Shown below) 
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Discussion 


Structural transition causing a conformational change 
occurs as a result of ligand binding, known as allosteric 
transition [9]. Antagonist such as THC, CBD, and propofol 
bind to different sites within the transmembrane domain. 
Binding of antagonists can create necessary driving forces 
for conformational transition in such a way that the 
closed channel state becomes more favorable. A strong 
binding affinity of antagonists causes more driving force 
for the receptor to stay in the closed state. Results of the 
lower binding site docking corresponds to _ the 
experimental finding that propofol is the strongest 
inhibitor of 5-HT3 receptor membrane currents (~50%) 
while expressing the largest binding energy of -6.0 
kcal/mol. The trend of absolute values of calculated 
binding energies for the other compounds 
(propofol>menthol> eucalyptol>cyclohexane) provide 
rational to experimental findings of difference in 
inhibition of 5HT3 by the same compounds [10]. Both 
THC and CBD strongly and selectively bind to the upper 
site (TYR 346). Similar binding energies for THC and CBD 
explain their similar inhibition of membrane currents in 
5HT3 [11]. In order to verify the predicted binding mode, 
we proposed experimental studies of site-directed 
mutagenesis of the residue involved in both hydrophilic 
and hydrophobic interactions at the THR 361 binding site. 
Experiments using thymol (-6.4 kcal/mol) and 2- 
isopropylphenol (-6.1 kcal/mol) will be done after 
docking simulations showed that both interact stronger 
than propofol and menthol with the THR 361 binding site. 
Thymol, a monoterpenoid similar to menthol, has an 
aromatic group which is responsible for its higher affinity 
to the lower binding site (THR361). 


l 00 7 (5 menthol alone 
* [ES compound alone (n=6) 
fon) HB mentho! + compound 
a =5 
Oo 80 =) 
—_ 
= 
— 
© 60} 
€ 
“= 404 
S 
20+ 


propofol 


cyclohexane —_ eucalyptol 
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Evaluation of the data in the figure on the left provided 
evidence that combining menthol with either eucalyptol 
or propofol enhanced its inhibitory effect on SHT-induced 
currents, suggesting that these compounds have additive 
effects on 5-HT3 receptors [4]. On the other hand, 
cyclohexane (100 mM) alone displayed undetectable 
efficacy at inhibiting 5-HT3A receptors [12], but increased 
significantly the potency of menthol, suggesting that 
cyclohexane may facilitate the interaction of menthol with 
its binding site on the receptor. 


Our computational findings did not back up the 
observed increase in inhibition of menthol when in the 
presence of cyclohexane, as shown in the figure above. 
Experimental error occurred while simultaneously 
docking two ligands in the lower active site (THR 361) to 
recreate these results. The initial process of docking 
menthol into the active site, saving it as a functional 
(flexible) part of the receptor, and docking the second 
ligand, cyclohexane, to the system of ligand and receptor, 
was proven ineffective. Cyclohexane, in the presence of 
menthol, causes a vast increase in inhibition than the 
compounds alone. This led to the theory that cyclohexane 
doesn’t pass into the transmembrane region, but sits on 
the outside of the channel affecting its membrane 
potential. The increase shown in the figure above could be 
attributed to the change in the potential energy across the 
membrane that cyclohexane could cause. 


Conclusion 


Experimental inhibition of homomerically expressed 
human 5-HT3 receptor shows similar trend in 
computational binding energies to the binding site located 
near THR 361. Docking calculations provide explanation 
of molecular basis of difference in inhibition by menthol 
like compounds. Similar binding energies for THC and 
CBD corresponds to their similar inhibition of membrane 
currents measured in experiment [3].These results 
provide rational for search of new compounds that might 
be used as antiemetic drugs. 
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